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ABSTRACT 


Under conditions of permanent gases diffusing through, 
foam bubble nenbranes with time dependent membrane surface 
area, a theoretical model for estimating gas permeabilities 
has been developed. For generating a monolayer of foam on 
the surface of a column of a liquid containing surfactants, 
an experimental glass apparatus has been fabricated. Experi- 
ments were conducted in this foam membrane perneator for 0 2 ~N 2 
and C0 2 -C 2 systems to determine the permeabilities of individual 
gases as well as the separation factors a q and . 

The relevant quantities like average bubble diameter, bubble 
shrinkage with time due to unequal permeabilities etc., needed 
in the theoretical model to find out the permeability of 
the gases have also been reasured. For C> 2 -K 2 s y s ^ em > three 
different compositions of the system glycerol - saponin- 
water were used by varying the glycerol content only and 70 
per cent glycerol solution was found to generate the most 
stable foam* 

The separation factors for the 0 2 ~1 2 system were in 
the range of 1.7 to 5.8 with 70 per cent glycerol solution 
having separation factors always greater than 3.5 because of 
improved foam stability. Those values are quite high and foam 
membrane separation of air is a distinct possibility. The 
separation factors for the CC 2 -0 2 system were in the range of 
6.7 to 8.3. These are somewhat lower than reported values 
due to possible absorption of C0 2 . 



CHAPTER 1 


i 

INTRODUCTION 

4 ^ ‘ ; 

Industrial separation of gases using membranes have 
recently generated considerable interest. Current applications 
invariably make use of a solid membrane along which the gas 
mixture to be separated flows at a high pressure and the 
permeating mixture of gases' having a different composition 
is withdrawn downstream of the membrane at a low pressure. 

Solid membranes have several disadvantages,- They have low gas 
permeabilities, and they are quite likely to have a few pin 
holes^ through which the feed gas mixture leaks to the permeate 
side without separation. In addition, capital investments 
are considerable while the operating costs are also high 
because of the high pressure differential needed to be main- 
tained across the membrane thickness. A possible alternative 
to solid membranes is provided by liquid membranes which were 
developed by Li [1] in 1968 for the separation of liquid 
hydrocarbon mixtures. The liquid membrane in this case is 
essentially a thin aqueous layer containing surfactants, which 
encloses a drop of liquid hydrocarbon feed and which in turn 
is surrounded by another hydrocarbon medium of large extent 
into which some species from the encapsulated drop will 
preferentially permeate through the aqueous liquid membranes. 
There are several advantages of such liquid membranes, whether 
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aqueous or nonaqueous (in the latter case the droplet and 
the receiving medium are aqueous). Their thicknesses can he 
very small and since diffusivities in liquids are in general 
much higher than that in solids, the permeabilities are 
likely to be much higher if the gases are reasonably soluble. 
Further the separation factors are likely to be much higher 
than that in solid membranes because of wide differences in 
solubilities of various species in liquids [2]. Further liquid 
membranes do not have pin holes. On the other hand one cannot 
have a high pressure difference across the liquid membrane; 

In fact, the biggest single obstacle to the use of liquid 
membranes for gas separation is how to support the liquid 
membrane and provide a large surface area necessary for 
industrial level gas separations. Thus even though liquid 
membranes are about to become commercially successful [3], for 
separation of liquid mixtures because of the simpler techniques 
of emulsifying droplets of one liquid mixture in another 
receiving liquid with an intervening, immiscible liquid membrane 
such is not the case for separation of gaseous mixtures. 

The earliest and the only reference to the use of a 
liquid membranes for the separation of a gas mixture is 
contained in the patent taken out by Li [4] in 1973. Li [4] 
separated mixtures of 1-G 4 H Q and 2-C 4 H Q by making a foam from 
a aqueous solution containing 70 per cent glycerol and 0.2 

The foam was carried up a column (of the 


per cent saponin. 
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"type normally used in foam fractionation) by a carrier gas 
(in this case l^). As the foam rises up the column,, driven 
by the surrounding carrier gas, the more permeable species 
will permeate at a higher rate through the aqueous surfactant 
layer into the carrier gas from which it can be separated 
later by absorption or liquefaction. It is obvious that in 
order to give sufficient contact time between the foam and the 
carrier gas the column has to be long and foam breakage becomes 
quite likely affecting the quality of separation. Further 
the carrier gas' velocity cannot be very high so that output 
rate is quite low. 

The motivating force behind this work is, to develop 
a foam membrane system for obtaining enriched air. Didwania[5] 
has already demonstrated the possibilities of an alternate 
way of using foam membrane system for gas separations by using 
foam in the form of a sheet carried down by a carrier liquid 
while exchanging gas with another carrier gas on the other side. 
In order to know what kind of separation is expected in such 
systems it was decided as a preliminary step to determine the 
permeabilities of these gases in a static apparatus using a 
foam membrane. Bulk of the investigations are intended for 
0 -* system which will provide more information for an air 
separation technique. The notable earlier investigation made 
in measurements of permeability of permanent gases, through 
flat surfactant films is that of R.L. Cook and R. W. Took 
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who used an aqueous film of 2 weight per cent Ivory liquid and 
2 weight per cent Dupanol for the formation of liquid film. 

The driving force for enrichment was achieved by maintaining 
a temperature gradient across the aqueous surfactant film 
since the concept of foam seemed to have escaped their attention. 

In this preliminary investigation an approximate mono- 
layer of foam is generated at the top of surface of a volume 
of liquid kept in a cylindrical glass container by bubbling 
a known volume of a gas species A through a sparger kept at 
the bottom of the liquid. The top of the liquid as well as the 
foam monolayer in the closed cylindrical vessel contains a 
second gas species B. The liquid membrane constituting the 
surfaces of the bubbles separates this second gas from the gas 
enclosed inside the bubbles constituting the foam. As time 
progresses, diffusion of the gas species B into the foam bubbles 
and diffusion of species A out of the foam bubbles take place. 

If the permeability of species A is more than that of the 
species B, the foam bubbles will shrink. One can simultaneously 
measure the decrease of height of individual foam bubbles 
with tine as well as determine the composition of the gases 
in the enclosed gas space on top of the foam bubbles. Unlike 
the situation in dialysis with solid membranes whose surface 
area remains constant with time, this case is one of membranes 
surface area and shape varying with time. The method of 
determining permeability of individual gas species is therefore 



5 


quite different and difficult from that in dialysis for 
solute species dissolved in two liquid streams separated by 
a solid membrane. The objectives of this thesis are therefore 
primarily as follows: 

1. Development of an appropriate analytical expression 
for predicting permeabilities from minimum measure- 
ments of concentration change and foam bubble height 
variation with time. 

2. Make preliminary measurements of the required quanti- 
ties in the following type of system: gas species A 
inside foam bubbles and gas species B outside with 
object of determining permeabilities of species A and 
species B and therefore the separation factor for 
A-B system with particular reference to the system 

W 

3. Obtain suitable information about foam stability as 
a function of the aqueous solution composition as 
well as about the nature of foam formed, the size 
distribution of foam bubbles etc. These informations 
are expected to provide an adequate guideline for 
desig nin g and operating a continuous foam membrane 
separator for a gaseous mixture. 



CHAPTER 2 


MODEL FOR EST IMATING GAS PERMEABILITY THROUGH FOAM! MEMBRANES 

In this chapter, a theoretical analysis of the bubble 
shrinkage due to unequal permeabilities in a static foam 
membrane, has been carried out and an appropriate analytical 
expression, has been developed for -predicting gas permeabilities 
from measurable quantities like gas composition change, foam 
bubble height variation with time, average bubble size etc. 

As pointed out earlier, separation of gaseous mixtures 
through foarn membrane is similar to dialysis in many respects. 

It is a rate governed membrane separation process. The 
process is passive without any external forces like pressure 
gradient or electrical potential gradient. No external 
energy is required for the separation of gases. The driving 
force for the diffusion of the gases is the concentration 
gradiont maintained across the liquid foam membranes. The 
problem of gas permeability determination with such an arrange- 
ment is, however, much more complicated than dialysis, since 
the surface area of the membranes is continuously changing with 
time. 

The system as shown in Figure 1, consists of a small 
column with liquid upto the level H^from the bottom with the 
space above being filled with a gas of volume Y . At the surface 
of the liquid a monolayer of foam is generated by bubbling 
another gas, whose permeability through the liquid membrane 








is to be determined along with that of the gas outside the 
foam. Initially the gas inside (gas A) and the gas outside 
the foam (gas B) are pure. As time proceeds, both gases 
permeate at different rates through the foam membranes it 
being assumed that the diffusional loss of species A into the 
liquid on which the foam stands is very little. This can be 
so only when the gases are sparingly soluble in the liquid. 
The species transfer process with time, leading to bubble 
shrinkage (if the gas A has higher permeability), can be 
described by the following equations: 

From a material balance of gas species A and B, 


dN A0 

at = “ 

N 

e 

3=1 

dn Aij 

dt 

(2.0A) 

dN_ _ 

BC _ 

dt 

N 

e 

4 d B15 

dt 

(2. OB) 


where 


3=1 


n, 




n 


Bij 


N 


AG 


N 


BC 


= moles ol species A inside jth bubble 
= moles of species B inside jth bubble 
s moles of A outside the bubbles (continuous phase) 
ss moles of B in continuous phase 


and N = total number of bubbles in the monolayer. 


If the moles of species A and B inside all the 
bubbles are denoted by and respectively, the above 
equations become, 
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and 


an 


BC 


dN 

Bi 

dt 


TT = " "dt” (2,2) 

Since permeation of various gas species is assumed 

to be strictly due to concentration gradient, the mass flux 

through tho liquid membrane around gas can be expressed in 

the form of Pick's first law, so that 

N 

e 

3=1 
N 

e S. Perm B[ 

*1 

3=1 


dK 


Ai 


N, 


n. 


dt 


‘"AC 1 ‘Aij 

S. Perm A [ P - — 2 -mJ 
3 \c C “Tij 0 


(2.3) 


dN. 


'Bi 


dt 


!bo p 

- L JJ £ C 

j "qjc 


_g.il p. 1 

n Ti j ^ 


(2.4) 


where P = total pressure outside the membrane in the 
0 

continuous phase 

p .= pressure inside the jth bubble 

1 J 

S. = surface area of the jth Bubble 

3 


N 


PC 


total number of moles of gases in the continuous 


phase 


n, = total number of moles of all gases inside the jth 
bubble 

Pern A *= permeability of species A 
and Perm B = pezmoability of species B 

To solve the above equations (2.1), (2.2), (2.3) and 
(2.4), the following assumptions have been made: 

(i) The gases are assumed to be ideal. 

v i,M n hrpakase (to prevent contamination) 
(ii) There is no bubble breakage vto p 
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(ill) There is very little absorption of either of 

the gases in the liquid phase on which the foam 
stands. 

Prom the ideal gas law, 

P, 

and 


■13 


n Tij 

V j 


(2.5) 


( 2 . 6 ) 


N T0 

p c = r~ 

(V- e v ) 

3*1 3 

where R = gas constant and 
T a* temperature 

and v^ = volume of tho jth bubble. 

The excess pressure inside tho bubble due to surface 
tension can be written an, 


2 cr 


*13 c 

whore a = surface tension 


r . 
3 


(2.7) 


r^ = radius of the jth bubble 

Since this excess pressure is very small (see Appendix 
A, 2.1) it can be neglected in comparison to P c . Therefore 
equation (2.7) becomes 



If Sj denotes the total surface area of all the bubbles 
through which species A and B simultaneously move in opposite 
directions, equations (2.3) and (2.4) may be written as 
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C ^Ai ^AC , . ^ ^Aij 

"dt~ * r~ p c Perm L S T^) - p c Perm A e S.(t) 

TC c . 1 3 n Tij 

(2.9) 


and dN B > N-gn ^ n -n 

-dt = P c B Sj^) - P c Perm B e Sj M ^ 

3=1 


n Bi.i 


F rom (2.5), (2.6) and (2.8) 


n Iij 


AT RT 

N 

P o (T- «T ) 


( 2 . 10 ) 


Substituting tliosc in (2,9) and (2.10), we have. 


dK [(RT Perm A) S T ] I AC B 

-# = B ~ - M Pcra A , e 

( V - c V.) 3 ~ 

3*1 °I IT 


N S.n, . . 
3 ^-3 


v. 

3=1 3 


( 2 . 11 ) 


[ (RT Perm B) S,] N BC S. n B±j 

a s RT Peru B e — 

( V - e v.) 3=1 J 


If wo write N Bi as Nn A± (t) and R- Bi = N n Bi (t) 
N 

e v. * N v(t), e S. = NS (t) 


( 2 . 12 ) 


S 1 n Bl-1 Cf h fLlUSA. (even though this is not strictly 

▼i " v 

J true but is likely to be true 

for moderately dispersed size 
distribution) . 
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where = moles of A inside the average bubble 
ng^ = moles of B inside the average bubble 

v = volume of an average bubble, 

S = surface area of an average bubble. 
Equations (2.1) and (2.2) can be written as, 


dN , „ 
il C 

dt~ 


dn 


N 


a i 


dt 


dN 


BG 


d 


N 


^Bi 


dt “ " dt' 

Addition of (2.13) and (2.14) will yield. 


(2.13) 

(2.14) 


dN 


IS. 


H 


d5 ri 


"dt™ = " " dt 

where n^ = + n Bi = total number of moles inside the 

average bubble. 


(2.15) 


Integrating (2.13) 


N - + K ss - N n, . (t) where K is constant of integration 

At t *s 0 , N n Ai (t) = N Ai which is equal to total moles of A 

inside all the bubbles at tine t=0. 

(2.16) 

V ' 'M' ■ 

Similarly, wo can got, 


Therefore (N^-N^) = - Nn Ai (t) 


(n bc - l^ c ) *= " N 
(»S - N tc ) - N ^(t) 


(2.17) 

(2.18) 


whore « total moles of B in continuous phase at time 


t » 0 
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and N 


T = total moles of A a nd B in the whole system at t=0» 


For equation (2.1l), 


Ratio I /XI = 


(RT Perm. A ) Srp N AC . 

~ - - - 

(V- e v )(RT Perm A) RS iL./v 
j=l 3 Al 


AG r } 

* if — — r - ^ ( Substituting for n, . 

(T- iv.) 

3=1 J from (2.16)) 

Here * i0 /<*£ i“^Ac) 1 for a very ; large time from t = 0 

(see Appendix A, 2* 2). 

Further if V»Mv , the term II in eqn. .( 2*11) is controlling. 
Hence equation ( 2. ll) may ''he simplified to 


dN . a 

-jr* « -RT Perm A e 

j=l 


S, n, 


(2.19) 


Substituting for -Nn^ ( t ) from-. (2.16) 
dN, iC (t) rt Pom A rt p , T , + n 


— siS . R? [K°.-R (t)] (2.20) 

dt V(t)/S(t) L ^ A0 

Solution of this equation for the initial condition t=0, 
N AC (t) = 0 and t=t , N IC (t)=known,. will enable us to determine 

Perm A. In practice, we have used known values of v(t), S(t), 
N AC (t) and [d ^(t )/dt] at any" given time to determine Perm A. 
The procedure for determining v(t) and S(t ) from r and h(t) 
aro given later in this section. The derivative is estimated 

from n ac (t) vs t data. 

Since the gas is assumed to be behaving ideally, 

N Aa . r m bc 

P AC “ P o Hjq ’ “ o H IC 
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n, . . 

p m , = p, , . p . = « 

xj ^Bij p 


n Bi.i 


U "ID n T .. 


where p stands for partial pressure. 

Substituting these in equations (2.3) and (2.4) 


and 


...N ...... A 

a5 iU 

N 

Perm A 

dt “ 

e 



3=1 

If 

dS Bi 

N 

Perm B 

dt = 

e 

3=1 


E S i ( PAC - PAij) 


( 2 . 21 ) 


( 2 . 22 ) 


Adding (2,21) and (2.22), we get 
d n. . w dH Bi 


N 


N ^ jfii jj 

Porn A TT" + Perm B ~d£~ e S j (p AC +p BG~ (p Aij +p Bij ^ ) 

3=1 

N 

Mi s i ( W 

which on using relation (2. 8) leads to 

dn^i Pern A dn^ 

"eft ” "* Perm B dt 

From equations (2.5), (2.6) and (2.8) we get 


(2.23) 


n Tii (t) N ic 
1 v'JtT " ~n 


d 


3=1 


v,(t); 


Assuming n Tii^i _ ^Ti^.l 




TC 


v 3 (t) 


v(t) 


N 

7- e v 
j=l 3 


we obtain using (2.18), 

Ajh (t) N **Ti 

wiJML mm 1 A 


v(t ) 


N 

7 - E v 
0=1 


3 



15 


*0?!^ 
v(t ) 


N 


Nn Ti (t) 


(2.24) 


Y- Nv(t) Y- Nv(t ) 

I II 

The ratio of two toms on the right hand side of (2.24) is, 

tO 


II 


N 


1 


Nn Ti (t) 


if wo assumed that V» Nv as we have already done before 
obtaining ( 2.19) j Nn^^(t ) so that the term II can be 

neglected in comparison to the term I. Therefore, equation 

(2.24) reduces to A 

H?V (t) 

JL r ♦ - 


n T 1 (t) 


j u 


dt 


dt 


V 

IJO 

"T 

~* L 

v(t) 


V 

expression 


dv 


[ * / v » Nv] 


(2.25) 


Thus (t) + n^(t ) 

Diff orontinting this expression with respect to t, we get 


( 2 . 26 ) 


An, 


UU* * 

Substituting for — from (2.25) one obtains 


ln Bj 

dt 


w° ‘_ 

r Perm B « JL. dv 
’•Perm B -Porn A-* Y dt 


(2.27) 


Integration with respect to time^leads to 

r f+.'s _ r -JlSES-S -1 v (t) + Constant. 

n^ti - Lp erm B-Perm A J V 


At tssO , n Bi (t) a 0 and v = v(0). 


Constant 


r Pem._3 j 

Lperm B~Perm A J 


Hi 


Y~ ^(o) 




^Bi* 


i+\ r ? erm B — -i -~ —f v(t ) -v( 0 ) ] ( 2 . 28 ) 

(t) * Lperm B-Perm A J * V L 
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From (2.12) 


N 

. ln Bi 

Porn B 

dt 



dn. 

<®»* «L 

Perm B 

dt ' “ 

; n bc ■ 

Q 

I 

H* 

1, 

dHj 

Pom 

B dt * 


K, 


R - C N s(t) - 


S(t) 

v (t) 


■] (2.29) 


F 


V f- N BC 


IZkL *2 

v( t) 


(since V» e v.) 

3=1 0 


■] 


S(t) RT 

V Li, BG ~ x, Bi 


o V^.(t) 

[tf n - ] (2.30) 


i * 0 . f 


v (t) 
pince R Bi = Fn^] 


clEr, 4 RT Perm B v 

_ r S(t) - -2- 


T 
dH Bi 


V (t) 


>] 


Substituting for and n B± (t) from (2.27) and (2.28) into 

this 

w° 

r norm B T (lv_ RT Perm B g( t ) ri P _r]^ -JL ] x 

ipcTS VptTnnTJ IT at ~ v sltJL ^ 1 z (t) 


£ 


ipsPPfe) 4- <**>-«»> 


dv 

i 

i • o . 


i: 


(R + r ) 


(2.31) 


BC 


II* s RT S(t) [(Perm B-PemA) v 


— — Perm B(v(t )-v(0 ) ) 


(F+ Z ) 


civ 

dt 


RT 13 (t) [(perm B-Penn A)*g Y Perm B(v(t )-v(0 ) ) ] 

(2.32) 


SX 



where x 


■B 


dv 


3C 


i* e *> dt" = ^ H (^ erm B-Perm A)x£ - i— Perm B (v(t)-v(O))] 

v(t) 


i . e . d v 
dt 


( See Appendix 2 . 23 ) 

RT S(t)[(perm B-permA)xg- perm B 

v(t) 

(2.33) 


The above equation can be solved to obtain perm B, 
knowing the initial conditions, the surface area and the volume 
of the bubble at any time t and perm A. This ordinary non- 
linear differential equation has been numerically solved after 
obtaining v(t) and S(t) as a function of r and h(t). The 
relationship between the surface area of an average bubble 
and its height h(t) as well as the radius of the bubble are 
dealt with next. 

For calculating the volume and the curved surface area 
of the bubble, its shape must be known. The bubble shape as shown 
in Figure 2, is assumed to be that of an oblate spheroid for 
which the volume and surface area can be calculated [ 6 ] as 
follows: 

The surface area of the average bubble is 

h 2 (t)r rt^ 2 -h^i 

s(t) = x ? 2 + ^ l 0 S e ~ 


V"(r 2 -h 2 ) 


v- 


(2.34) 


where r - radius of the average bubble which is assumed to be 


constant. 
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h = height of the average bubble, which varies as 
„ "the gases diffuse through the bubble 
The volume of the average bubble is given by 

v(t) = j % r 2 h(t ) (2.35) 

Hence by measuring the radius and height of the 
bubble wo can calculate the surface area and volume of the 
bubble at any time t. These values of the surface area and 
the bubble volume can be used to estimate (dv/dt) for use in 
the equation (2.33) for determining permeability of species B, 
since Perm A, B(t), v(0) and v(t) are known. 

From the permeabilities of A and B which can be obtained 
from the equations (2.20) and (2. 33,)* respectively, we can 
calculate the ideal separation factor which may be defined as 


Separation factor = 


Perm A 
Perm B 


(2.36) 


An absolutely rigorous procedure would have been to solve 
the nonlinear ordinary differential equation (2.33) for h(t) 
and also obtain an analytical solution for N AC (t) from (2.20). 
Since h(t) and N AC (t) are available experimentally, Perm A 
and Perm B could have been estimated. 


The permeabilities of the gases which can be determined 
from tho equations (2.20) and (2.33) are largely governed by 
the solubility of the gases in the foaming liquid. But it can 
bo improved by adding a facilitating agent which enhances the 
permeation of one of the species [If ]. The permeability values 


for the facilitated membrane has bec.n given in Chapter 4. 





CHAPTER 3 


EXPERIMMTAL DETAILS 

This chapter is divided into three sections. The 
first section deals with the chemicals used for the preparation 
of the foam membrane. The second section deals with the 
description of the experimental set up for gas permeability 
measurements and the third section explains the procedure 
used for taking data for the calculation of the permeability 
of gases A and B and the separation factor g ^ . 

The steps followed during the course of experiment 
were as follows: 

(i) Preparation of foaming liquid and filling the 
liquid in the glass column. 

(ii) Formation of the foam membranes by bubbling gases 
through the sparger. 

(iii) Measurement of the various factors used for the 

■ . -*v 

permeability calculation. 

(iv) Breaking the foam membranes by draining the liquid 
and cleaning the apparatus with water. 

3.1 Chemioals : 

The chemicals needed for the preparation of the 
foaming liquid are glycerol, saponin, haemoglobin and water. 
Saponin acts as a surfactant for generating the foam. The 
glycerol used was a laboratory reagent grade while saponin 
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was of imported analytical reagent grade. The other chemicals 
needed were sodium dithionate and KOH for oxygen content 
measurement. Sodium sulphate, and methyl orange were 

used for the solution in the aspirator hottle of the Orsat gas 
analysis apparatus. For the absorption of oxygen 80 gms of 
sodium dithionate was dissolved in 40 ml of water. This was 
then mixed with 48 gms of KOH dissolved in 40 ml of water. 

The solution was shaken thoroughly before using for the absor- 
ption of 0 2 in the Orsat apparatus. The solution used in the 
aspirator bottle consisted of 20 gms of aqueous sodium sulphate 
dissolved in 80 ml of H^O and, 4 gms of H^O^. Methyl orange 
was used to give colour to the liquid. 

3.2 ' Foan Membrane Permeator Details : 

The design basis of the apparatus is available in 
Appendix A. 3.1 at the end of this thesis. The apparatus (Figure 
l) consisted of a corning glass-chamber of 6 inch diameter and 
20 inch length. The chamber had a glass sparger of 5 inch 
diameter with 50 orifices of 1 mm diameter. The size of the 
orifices in the sparger had to be uniform to get bubbles of 
uniform dimensions. Through the stopcock (l) , the chamber 
was connected to the gas cylinder containing gas B, while 
through the stopcock (2), the chamber was connected to the 
manometer filled with mercury. The gas to be analysed in the 
Orsat apparatus was taken through stopcock (3), which is a 
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two-way high vacuum stopcock. The foaming liquid kept in 
the feed bottle at a height to provide sufficient head 
for the flow of the liquid, was sent through the stopcock (4) 
while the draining had already been made through stopcock (5). 
The sparger kept inside the chamber was connected to a high 
vacuum three-way stopcock (6), the other two ways being 
connected to the vent line and the wet-test meter, which was 
connected to the gas cylinder containing gas A. The entire 
flow loop is shown in Figure 3. 

3 . 3 Experimental Procedur e: 

Tho foaming liquid was prepared with various concen- 
trations of glycerol to find out the effect of liquid compo- 
sition on separation factor. It is assumed that increased 
glycerol concentration will make the foam stabler since Lite [4] 
liquid membrane stability increased very much with increase 
in glycerol content. For 0 2 «N s 7 s "k en > concentration 

the foaming liquid used were: 

(i) 30 per cent glycerol, 68 per cent H^O and 2 per 
cent saponin. 

(ii) 50 per cent glycerol, 48 per cent H 2 0 and 2 per 
cent saponin. 

(iii) 70 per cent glycerol, 28 per cent H^O and 2 per cent 
saponin. (without facilitating agent) 

(iv) 70 per cent glycerol, 28 per cent H 2 0 and 2 per cent 
saponin (with Haemoglobin acting as a facilitating 
agent ) . 
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Per CO 2 -0 g system the liquid used was 

70 per cent glycerol, 28 per cent H 2 0 and 2 per 
cent saponin. 

Foaming liquid thus prepared was fed into the chanter 
from the feed bottle shown in Figure 3. Sufficient space 
was provided at the top of the foaming liquid in accordance 
with the assumption Y"^>nv, made in the development of 
theoretical equations in Chapter 2. The space above the 
liquid was filled with gas B, sent through stopcock (l) , 
after flushing the space with gas B. Minimum pressure (about 
1.01 atns) had been kept at this space to minimise breakage 
of bubble during the bubble forming process. The system was 
ready for operation under such conditions. A monolayer of the 
foam of gas A was formed on top of the foaming liquid, by 
slowly sending the gas A through the sparger which in turn 
was connected to tho wet-test meter. Because of the higher 
permeability of gas A, the foam bubble started shrinking with 
time immediately. The various quantities needed for the 
calculation of permeabilities of gas A and B (and hence 
separation factor) are: (a) total moles of gas B initially 
present, (b) total moles of gas A supplied for making foam 
(assuming no loss); (c) the average radius of the bubble, 

(d) height of any shrinking foam bubble as a function of 
time; (e) radius of the bubble for which the reduction in 
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height was measured; (f) composition of the gases in the 
continuous phase at various times, as it changed due to the 
pemeation of A and B through the foam membrane. 

Total moles of gas A at tine t=0 was obtained by noting 
the pressure in the continuous phase and the volume of the gas 
space above the foaming liquid. The pressure was noted with 
the help of a nanometer attached to stopcock (2), while the 
volume was calculated from the height of the gas space and 
the radius of the chamber which were measured. The total 
molos of gas B at tine t=0 was obtained from the pressure and 
the volume of the gas fed in through the sparger. The volume 
of the gas A fed in was obtained from the wet-test meter, while 
the pressure was noted from the manometer connected to it. 

The average radius of the bubbles was obtained by photographing 
bubbles formed on the top of tho foaming liquid. The total 
number of bubbles in a particular area of the photograph was 
counted to get the size of the average bubble. The height 
of the bubble at various interval of tine was measured by 
focussing the bubble to the crosswire intersection of a catheto 
meter. Every minute the crosswire intersection was made to 
coincide with the tip of the bubble and the corresponding 
reading was taken from the scale attached to the cathetometer. 
The radius of the bubble for which the height reduction was 
noted was obtained from the angle of the arc and the radius 
of the arc of the circle to which the diameter of the bubble 
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was the chord. The angle was measured with the scale provided 
at the bottom of the cathetoneter. The radius of the arc, 
which was the distance of the bubble to its image in the eye- 
piece, was obtained by focussing an object without disturbing 
the eyepiece and measuring its distance from the image. From 
the angle and the radius of the arc, the chord length was 
calculated. This was the diameter of the bubble. Composition 
of the gases in the continuous phase was measured by taking 
75 cc of the gas from the continuous phase and analysing it 
with the help of Orsat apparatus. From the reduction in volume 
of the gases in the gas measuring burette of the Orsat apparatus 
after absorption of one of the two species in the continuous 
phase, the percentage by volume of species A and B was calculated. 
Since volume per cent = mole per cent, these values were 
substituted in the equation (2.20) for the calculation permea- 
bility A. 



CHAPTER 4 


RESULTS ADD DISCUSSION 

.In this chapter, the results of the experiments 
carried out to determine the permeability of gases through 
foam liquid membranes have been presented. Experiments were 
conducted primarily for the system with various compo- 

sitions of the water-glycerol saponin foaming liquid. A few 
ru is were made with the CO 2”^ 2 system also. The system Og-Hg 
is important for separation problem and therefore extra 
emphasis was put on it to the exclusion of the other systems. 

Due to the nonavailability of large number of permanent gases 
and the difficulty in analysing the compositions of such gases, 
only limited data with limited number of other systems could 
be taken. The results of the experiment conducted for the 0 2 -D2 
system with 30 per cent glycerol solution are not reported here 
since the extent; of foam breakage was found -to be large because 
of higher foam drainage rate with low viscosity solutions. 

The results are presented under the following headings; 
(,i) Effect of liquid composition on foam stability 

(ii) The nature of foam bubbles 

(iii) The variation of foam bubble height with time. 

(iv) Determination of gas permeabilities. 

(i) Poam Stability; Visual observations of monolayer foams 
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produced in the experiments indicated that as the glycerine 
concentration increases the foam stability increases. Rough 
counting of the number of foam bubbles breaking in a given 
interval of time was done and as the glycerin content of the 
2 per cent saponin solution in water increased from 30 per cent 
by weight to 70 per cent by weight, bubble breakage became a 
rai’e event. Only three compositions were usea: 30 per cent, 

50 per cent and 70 per cent glycerol by weight. Experiments in 
the higher concentration of gLycerol could not be carried out 
since the handling of the solution posed a serious problem with 
highly stable foams being generated during any transfer operation 
between reservoirs. The stability was in general independent 
of the bubble sizes except in the case of very small bubbles 
having diameters around 0.20 cm. These were much more stable 
than the bigger ones with even 30 per cent glycerol solution. 

(ii) Nature of Foam Bubbles ; A schematic diagram depicting 
the actual nature of foam bubbles from a position right above 
them is shown in Figure 2. A photograph of the bubbles is also 
attached herewith (Figare 4A; . It may be noticed that the 
bubbles less than 2 cm in diameter but greater than 0.6 cm. The 
others had size ranging between 0.6 cm and 0.1 cm. However, 
their total surface area was almost always less than 5 per cent 
of the total surface area of the bubbles. No attempt was. made 
to study the height reduction of such tiny bubbles. Visual study 
of the bubbles further revealed that initially (t=0) most of 
the bubbles were of hemispherical shape. There is always some 
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FIGURE 4a: PHOTOGRAPH OF THE BUBBLES 



30 


coalecence at the beginning of each experiment. As time 
progressed, it changed to a form similar to an oblate spheroid 
with the height reduced. This supports the assumptions made in 
Chapter 2. 

3. Variation of Foam Bubble Height with Time ; In Figures 4,5,6 
and the 7 the variation with time of the height of several 
foam bubbles are shown. In anyone of these figures, the height 
variation with time for a particular bubble is very rapid in 
the initial period, but becomes almost zero at a large time. In 
so far as height variation with time is indicative of the 
separation process, obviously the rate of separation decreases 
with increasing time since the concentration driving force 
decreases with the concentration of A increasing outside the 
bubbles. It may be pointed out .also that even though anyone of 
the figures 4,5,6 and 7 contain a number of curves for different 
sized bubbles, only one curve can be generated in a given experi- 
ment due to the limitations of one observer and one cathetometer. 
However, experiments were repeated under conceivably identical 
conditions for any particular system and the height variation of 
a similar bubble was noted against time. Since the curves in a given 
Figure (e.g. Fig. 4 or Fig 5 or Fig. 6 ) are for bubbles with 
different size and the theory for this work as pointed out in 
Chapter 2 requires the height variation of an average bubble 
whose size is estimated from photographs, the height variation 
with time of such an average bubble has been obtained by 
Lagrangian interpolation formulae [9] (see Appendix A.4.1) and 
tabulated in Tables 4.1, 4.2, 4.3 and 4.4. 




FIG, 4 

BUBBLE HEIGHT VARIATION. WITH TIME FOR 02~N2 SYSTEM 
( 50% Glycerol ) 



FIG. 5 

BUBBLE height variation WITH TIME FOR 0 2 -N 2 SYSTEM 

<yd7o- Glycerol ) 
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HEIGHT VARIATION WITH TIME FOR C0 2 -0 2 SYSTEM(70*/. Glyce- 
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(iv) . Determination of Gas Permeabilities : 

The various data needed for the gas permeability 
calculations are tabulated in Tables 4.1, 4.2, 4.3 and 4.4. These 
quantities were substituted in equations (2.20; , (2.34; , (2.35) 
and (2.36) and solved numerically to determine Perm A and Perm B. 
A model calculation for determining gas permeabilities as well 
as the separation factor ' is shown in Appendix (A.4.2) . The 
results are tabulated in Tables 4.1, 4.2, 4.3 and 4.4. 

A look at Tablds 4.1, 4.2, and 4.3 shows that the sepa- 
ration factor for the oxygen (A) - Nitrogen (B) system 
increases general with (l) increase in gLycerol concentration 
a.nd (2) with time. Further the separation factor values are in 
the range 3.4 to 5.8 with 70 per cent glycerol solution. In 
addition, the introduction of haemoglobin, a known 02 - facili- 
tating agent, does not substantially increase the separation 
factor with is quite unusual. The explanation for this type 
of results are to be sought from the experimental observations, 
nature of data reduction and the numerical solution techniques 
adopted. 

dDhe addition of 0.3 per cent haemoglobin to the 70 per 
cent glycerol solution did not result in a true solutioneven 
through 0.3 per cent protein concentration is quite low. It was 
found that a large amount of "haemoglobin was floating on the 
surface of the foaming liquid. Thus most of it was not available 
for facilitation. This must be due to the large gLycerol 
concentration in the particular solution (70 per cent) in which 
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haemoglobin was used. Further it was found at the later stages 
that the haemoglobin used was denatured. However, the large 
time separation factor for this system is Quite appreciable 
(around 5 ) . . 


The increase in oc^g with increase in glye -rol concen- . 
tration in the foaming liquid is most likely due to the sharp 
reduction in bubble breakage with time resulting from a radically 
increased foam stability. However, there are other factors here. 
An increase in glycerol concentration will also lead to a thicker 
foam membrane which will reduce gas permeability (of all gas 
species) . Consider Table 4.1 and Table 4.2, the 0 ^ permeability 
as well as IT permeability decrease with increase in glycerol 
concentration from 50 per cent to and 70 per cent at low times. 
This is most likely attributable to thicker foam membranes. This 
phenomena was observed bjr Li [5] as well for liquid membranes 
separating liquid mixtures. 

The change in the values of the gas permeabilities as 
well as separation factor with time is to be looked at from 
several directions. The data taken at low times will be affected 
by initial bubble breakage that takes place during bubble 
formation. This reduces a ^ as well as permeabilities because 
of contamination. Such is not the situation for large time with 
stable bubbles. In addition, as time proceeds the foam drainage 
in the bubbles will considerably reduce the foam membrane 
thicknesses leading to higher values of permeabilities. Table 4.2 
and 4.3 demonstrate this. The increase in separation factor ct^g 
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with time can only be partially explained by lack of bubble 
breakage which takes place at the time of forming bubbles. An 
additional explanation is to be -sought from the manner in which 
the term (dN^/dt) was calculated. This term, was approximated 
by ( £. 3T^q/ At) and experimental values substituted to obtain 

its value. However, in the initial time period, the height 

■ , ■ .■>■■■■■ - 

variation becoming almost zero and any measurement error will 


bo considerably magnified. 

It is quite appropriate to compare the a Q values 

of this work with those available in literature for separation 
through flat liquid membranes. Oook and ‘lock [8] , he„d obtained 
a value of 2.24-2.99 for system with 2 weight per cent 

surfactant solutions in water. N.jff.ii [ll] had obtained a 
separation factor of 4 for Og-Hg separation through pure glycerol 
film. The present result of between 3.4 to 5*8 is 

therefore understandable in such a context. 

The results for CC^-Og system as shown in Table 4.4 are 
similar with regard to increase in permeability with time. However, 
the separation factor decreases with time. In addition, the 
magnitude of this separation factor a nr , is around 8, 2-6.7 

which is smaller than that reported by Cook and Tock [8] (namely 
19.9) . The most likely reason for this behaviour is that a lot 
of CO 2 is likely to have been absorbed by the foaming liquid which 
was not in equilibrium with COg. Further, CC >2 permeabilities 
being higher, measurements have to be done quickly. This had 
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introduced unexpected errors since only 3«5 minutes were 
available for the experiment. If CC> 2 were used, as gas B after 
equilibriating the foaming liquid with it and 0 2 were used to 
generate bubbles, some of these difficulties could have been 
avoided. ' ' 

It must be clear from the results quoted and the 
discussion that a significant advance in. gas separations is 
possible through foam' liquid membranes. , However,, the present 
measurements of gas permeabilities through foam liquid membranes 
have the following inaccuracies: 

1. lack of simultaneous measurement of height reduction 
with time of a number of bubbles in any experiment, 

2. Orsat gas analysis methods ar§; slow and chromatographic 
methods are necessary., 

3. Sparger " design should enable the generation of a 
monodisperse bubble population, 

4. The foam membrane permeator column surface was 

improper resulting in difficulties for obtaining 
a sharp photograph. 



CHAPTER 5 


CONCLUSIONS 

1. It is possible to make reasonable estimates of 
permeabilities of permanent gases through liquid membranes 
in the form of a foam liquid membrane and thereby determine 
the separation factors for separation of gas mixtures. 

2. Foam, stability is quite important in making such 
measurements accurately. 

3. The selection of the gas system and the liquid for 
foam making must take into account gas absorption or prior 
gas-liquid equilibration. 

4. Very high separation factors of around 5 are possible 
for the Og— Ng system with 70 per cent glycerol solutions whereas 
solid membranes have Og-ITg only in the range of 2-3. Thus 
continuous foam membrane separators are likely to become 
possible candidates for air separation as well as other 
separations. 

5. Foam membrane gas separation processes with monolayer 
foam in the sheet form has been shown to be an altogether new 
approach in gas separations. 



CHAPTER 6 


SUGG-ESTIOUS' FOR FTTfURS WORTT 

’ i 

X. More accurate methods for determining the composition 
of the gases in the continuous phase have to be devised. This 

ay be done by means of a thermal conductivity detector a 
Ga s-L i quid Chroma.t o graph . 

2. The surface of the glass column in which the foam is 

generated has to be smooth so that good photographs may be 
obtained. 

3. A permanent stand for keeping the camera with special 
close up lens, has to be made. 

4. Special lamps have to be used for uniform illumination 
of the bubbles to obtain a sharp photograph. 

5. Provision has to be made for breaking the bubbles, 
inside the permeator and hence cleaning the apparatus. This 
can be achieved by keeping a metallic ring inside the glass 
column, which can be moved from outside. 

6. The holes on the sparger should be made uniform and 
symmetrically to got uniform bubbles in a monolayer. 

7. The simultaneous measurement of height reduction of 
several bubbles can be made with a movie camera. 

8. Gas permeability measurements with A and B inter- 
changed should be at tempted, 

9. Attempts should be made to obtain analytical solution* 
for doe non-linear differential equations governing the 
systems. 
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10, With closer data points visa-a-vis time, better 
estimates of the first derivatives can be made, 

11, Explore the separation of a gas mixture inside the 
bubbles with a rather impermeable gas on the outside. 

12, The thesis can be extended to a continuous system 
as discussed by Didwania [5]. 
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APPENDIX A. 2. 2 


ESTIMATE 01' MOLE RATIO OP SPECIES A IN BOTH PHASE' AT 

LARGE TIME 


Approximate valu 
When the time is 


e of Kx = 


large. 


0.02 moles 


approximate value of = 0.011 (say) 

^Ai^AC = 0,02 “ °* 011 


0.009 





0.011 

0.009 







APPiSjffDlX A i 2 i 3 



H + - — 1 — : 

V(t) 

I II 

1 + II = II + 2 OH 


10-gjO) 

vT 6)/2 


Since the second tern is 20 times larger than th , s f±rst 
it can be approximated as N + - X a 7 V 

V(t) ?(tj 




APPEEDIX A.^-l 


BASIS 0P_. DESI G EIEG FOAM M8MEP ANE PERMEATOB 


Prom Pick’s law 


Q A + I K A P (-AA- *AB 


X, 


) 


(A. 3.1) 


Q 


'A Flux of th ~ ms A in (gm. moles/sec) 


v = Film thickness (cms) 
P = Pressure (Torrs) 


K 


L A 


"AA. 


— Permeation constant [-^P 1 * mo les -1 

Sq.cm. Torrs sec- 1 

Concentration of gas A on sido A 


C 


'AB' Concentration of gas A on sido B 
Knowing Q A , P, K^ 9 w and concentration of A both inside 
out old o, A can bo calculated. ¥e now consider a gas like 
C0 2 or 0 2 for estimating the surface,; area required. It has beer 
observed that 50 per cent of Go' diffuses out from the bubbles 
within 3 alii, from starting ( 5 ) with C0 2 inside and B- outside. 

Bow the volume of gas space above foaming liquid = 4330 c 
Volume of gas diffused = (4330)4/100 

= 173.2 cc 

. ,gn.aole S of gas diffused gm.moles 

q a = %kkm§ 

assume a film thickness (10) of 100 A 0 ; 

w = (100,10 

Pressure = 1 atm = 760 Torrs. 


.moles 
sec 


“ 8 cm = 10“ 6 cm. 
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Fron (8) , E. = (12.3) 10“ 9 - cc « m* noles 

sq. cm. Terraco 


Suppose we assume x = 0.2 and 


X AB = 0 


then 


^AA ~ 0*2(760) Torrs. 


Substituting in eqn. (A.3.I; 


... 173.2 (10 r\~9/ A \ 

22400^1807“ = ~ fb ° (A) * 760(0.2) 


A - - 173.2 

O /t 1 1 *ST" ; 


>7P0.2 10" 


= 0.2 x 10"* 4 cm 2 


Hence oven for 10~ 2 cm thick bubbles, surface area 

heeded is very little. For sake of convenience and to got 

.-epic gas fo. analysis wo make thts column of 65" dia. and 20" 

Height. This Will reduce, the edge effects on bubbles and 

fairly wide area would bo available. Fo r Orsat analysis at 

least 4 or 5 cc of p erne at or gas would b- ^ 

6 wouia Do available per 100 cc 

carrier gas' (Phase B) . 




APPEHDIX A, 4.1 


1, AG-RAIm'GIaN INTERPOLATION FOR E(t) 


A model calculation, of tho Lagrangian Interpolation 
necessary for estimating E(t) is shown below. 

lor 70 per cent glycerol, 28 per cent H^O and 2 per cent 
saponin solution, tho r and h values at time t=5 min were 
obtained as follows: 

r(cm) = 0.570 0.680 0.732 
h( cm) = 0.265. 0.393 0.513 

From photographed data r = 0.590 cm for 70 per cent glycerol 
system. The value of E(t) ^ minutes corrGS P on< ^ in g 'to this 
value of r is to be interpolated from Lagrangian interpolation. 


By Lagrangian interpolation formula 

^(■*•7 - Pi & 3 P 2 (A.4.l) 

The value P^,,P.^ and P^ may be found by 
(r-r ? ) (r-rj 
o ( r i** r 2 j (r-r) 


- ( 0 *596-0. 68) (0.596-0 .732) 
“ (0 .57 0-0 .68) (0 . 570^677327 

Similarly P^ and P^ may be written as 

P - (0*596-0.57) (0.596-0.732) 
1 ~ (0.680-0. 5a) (0.63-0.732) “ 


and 


p = (0.596-0.57) (0 596-0.68) 
2 ( 0 . 732-0. 5 1 HO . 732-0. 68) - 


0.641 

= 0.618 
-0.259 


These values when substituted in (A.4»l) 

E(0.596) . ’ = 0.265(0. 64l)+(0. 396) (0.618)-(0. 513) 

5mxnut e s 

(0.259) 

- 0,2727 



APPENDIX A. 4. 2 


MODEL CALCULATIONS POR DETERMINING PERMBA3ILI1Y OP SPECIES 

A AMD 3 


SYSTEM: 0 2 -F 2 .... 

EOAM LIQUID GOUGE ORATION: 50 per cent Glycerol, 48 per cent H 2 0 
and 2 per cent Saponin. 

Height of the column above bubble surface = 18.2 cm. 


Dimeter of the column 

Volume of the space filled with N 2 

Pressure in the continuous phase 


=17.0 cm 
= 17 2 x 18.2 x -re/4 
= 4131.0 cc 

= 1.022 atm 

(contains only N ? to 
start with) ^ 


Moles of F 2 in the continuos phase at t = 0 

= PV/RT 


Moles of U 2 , N° c 


R = 82.06 atm/gm. moles °K 

T = 308 °K 

_ 1.022 x 4131 
82.06 x 308 

= 0.167 gm. moles 


Volume of 0 ? sent in to form foam bubbles = 400 cc 


Pressure of 0 2 
Moles of 0 2 j 




= 1.03 atm 

_ 1.03 x 400 
82.06 x 308 

= 0.0163 gm. moles 

= WMMi 


= 0.911 

Radius of the average bubble from the photograph = 0.371 cns. 
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At tine t=0 nin, volume of the average bubble v(0) = h ? 2 z 2/3 

' ■ = 0.1062 cc 

At tine t - 3 min, volume of average bubble v(3) = 0.0797 cc 
Surface area of the average bubble s(t)= ; 

r 2 + iLtLJs lh r+(r 2 -E g )* . ... .. 

(r-h 2 ) 2 r 2 -(r 2 -E 5 )^ 

s (3) = 0.7224 cm 2 


If we assume the change in pressure and volume in the 

continuous phase is neslisibl e j /ii-iei 

xiegrigiDic due bo difiusion of species A 

and J3, the total moles remain the same. 

Concentration of 0 ,, at t-3 min = 6.5 per cent 

l'T A0 (3 min) = 0.065 x 0.167 = 0.0108 gm.moles 

= ^ Peria - A (^i-\ C (t) ) 


dP 

ac 

dt“ 


w(t;/s(t ; 

— 0.0036 gm. moles/min 
0.0797/0.7224 cm = 0.110 cm, 


iL 3) 

5 ( 3 ) 

“Aif’AoC* “ib) = 0.0163 - 0.0108 ='0/0055 gn.ooles 


Perm A 


dv 

dt 


_ 0.0036 x 0.07Q7 

0/7224 s 82.06 x 308 x 0.0055 

= 2,8 x 10" 6 , moles 

(atm; (cm* ) min 

= Pis (Perm B - Pern a) x^-Perm B( l-v(o) ) 


- 7 2 dE 2/3 

dt 

= * 0.371) 2 (0^-3685-0.2765) , . * 

m- . 0-3 (2/3) 

— —0.0088 cc/min. 


v(t) 



55 


v(0)/v(3) = 0.1062/0.0797 = 1.3341 

0.0088 = (82.06) (308) (0.7224) ((Per* B-Pe m A) 0.911 
Pern B (1-1. 3341)) 


i.e. , 1.2451 Pera 'B-o.911 p er n A = -0.48 x 10~ 7 

Porn A = 2.8 x 10 gn. nole S /(ci:i 2 ) (atn) (nin) . 
5-n B = 1. 66 x 10 -6 ^ nolea/m2 ‘ ?tc ; nto> 
S«p arator factor cc — — -A. _ _ ^ 


: AB ~ PerrTB” = 2 ‘ 8 * 10~ 6 /l .66 x 10~ 6 = 1.687. 




ITOME JTCL ATTIRE 


Tho following are moat commonly used symbols in the theis. 
n Aij = moles of species A inside jth bubble. 
n Bij = moles of species B inside jth bubble 
If A0 = moles of A outside the bubbles' (continuous phase) 
hsc “ moles of B in continuous phase 
H = total number of bubbles in the monolayer 
? 0 = total pressure outside the membrane in the continuous 


i j 


S 3 

S 


phase 

= pressure inside jth bubble 
= surface area of jth bubble 


*C = t0tal numb0r 0f “ 0lcs of *“»<> in the continuous phase 

= total number of moles of all gases inside the jth bubble 

Perm A- permeability of species A 

PermB= permeability of species B 

P- ~ gas constant 

P = temperature 

= volume of tho space above liquid level 
= volume of the jth bubble 
= surface tension 
= radius of the jth bubble 
= total surface area of all the bubbles 
n Ai = moles of A inside the average bubble 
n Bi = molGS of B inside the average bubble 
= volume of an average bubble 


V 

V j 

c 

r . 


S, 


■T 



S surface area of an average bubble 

U Ai = total moles of A insido the bubbles at time 
P. = partial pressure 

r = radius of the average bubble 

E(t) = height of the avorago hubblo at any time t 
a AB = se P ara -tion factor 
A = species A 
B = species B 


Subscripts 


i inside the bubbl s 

3 = jth bubble 

c = continuous phas^ 

T = total 

Supurscrl pts 

= average 

o - values at time t=0 



